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Abstract 

The formulation of the Lifshitz formula in terms of real frequencies is 
■ reconsidered for half-spaces described by the plasma model. It is shown 

^sO ! that besides the surface modes (for the TM polarization) , and the photonic 

modes, also waveguide modes must be considered. 



\Q ■ PACS: 12.20.-m, 12.20.Ds, 41.20.Jb, 42.50.-p, 73.20.Mf. 

1 Introduction 

In this paper the vacuum energy of the electromagnetic field in the presence of 
^ ■ two half-spaces characterized by a permittivity e and separated by a gap (see 

Fig. [Q, will be considered in terms of the frequencies uj of the physical modes, 
symbolically 

J 

This is the vacuum energy, defined as the half sum of the energies of all excita- 
tions, numbered by an index J, in the sense of zero point energy following the 
approach used by Casimir in pQ. At once, for the considered configuration, this 
is equivalent to the well known Lifshitz formula [2], which is usually written in 
terms of imaginary frequencies, u = i£. 

There are the following reasons to do this. 

First, the corresponding mode sum (see Eq. ([3]), below) was never written 
down correctly for the plasma model. The point is that, besides the modes 
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corresponding to surface waves and traveling waves (photonic modes), in ([I]) also 
waveguide modes must be included. 

A second reason comes from a situation where the permittivity e(u) is fre- 
quency dependent and known for real frequencies u only (e.g., from experimental 
data). In such case, the analytic continuation to e(i£) may introduce additional 
complications and a real frequency representation is of interest. 

A third reason comes from the discussion of the contribution of surface modes 
to the Casimir effect a few years ago in [3j 0] which might need to be reconsidered. 

The Lifshitz formula, written in terms of imaginary frequencies, has the signifi- 
cant advantage to involve only fast converging integrals. In contrast, integrations 
over the real frequency axis are notoriously difficult to handle because of the 
inherent oscillations. It will shown how this problem can be handled. 

In the following we consider a frequency dependent permittivity as given by 
the plasma model, 




*(«) = l-^> (2) 

where u p is the plasma frequency. For the considered configuration, the regu- 
larized, but yet unrenormalized (equipped with a tilde) expression for vacuum 
energy, following from (pQ), reads, 

- 2s + r ~^- 2s ^s}, (3) 

where s > | is the regularization parameter (in fact, Eq is the zeta function for the 
considered configuration). The modes include, for the TM polarization, the two 
surface modes with frequencies u s f, and, for both polarizations, the waveguide 
modes with frequencies Uj (j = 1, . . . , [— ] + 1) and the photonic modes with 

frequency tu > \J^l + &jj (in units with c = 1). The latter frequencies constitute, 

for fixed k||, the continuous part of the spectrum and enter with a weight given 
by the scattering phase shift S (its derivative is the mode density), see Eqs. (H3|) 
and (l44j) . The vector k|| denotes the momentum in parallel to the interface. The 
details will be given in Section 2. It is natural that the vacuum energy needs for 
a renormalization. It is only the force, 

f = ~Tl e " (4 > 

which is finite by itself (however, only after subtraction of the empty space contri- 
bution). In Section 2 we will see that the renormalization gives a useful suggestion 
for the treatment of the sums in ([3]). 

The vacuum energy, following after renormalization from ((3]), must coincide 
with the known representation in terms of integration over imaginary frequencies 
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which can be written in the form, 



_ 1 f dk\\ f°° di d ( 1 1 \ 

with the reflection coefficients, 

k-7] x-e{ii)i] 

r T E = — — , r TM = — ■ — 7TT7-, (6) 
x + 77 x + £(24)77 



and the notations x = y/e(i£)£ 2 + fey, 77 = y/£ 2 + A; 2 . The integrations in flHJ) 

converge fast due to the exponential factors. Therefore no renormalization or 
regularization is necessary in this formula. Also, it can be seen, that it is normal- 
ized in the correct way as to be vanishing for infinite separation. In fact, for large 
separation L, one can put the reflection coefficients (jHJ) equal to unity (formally 
by uj p 00) and will recover the result originally obtained by Casimir in pp. 

Eq. ()5]), or the force derived from it by means of (HI), is one of the many ways 
of representing the Lifshitz formula in terms of imaginary frequencies (see, e.g., 
Eq. (12.29) in [5], for reference). Other representations of this kind differ only by 
changed notations, by substitutions of the integration variables, by integrating 
by parts in £ or by considering finite temperature T. The latter can be achieved 
by substituting the integration over £ by the Matsubara sum. 

In general, the Lifshitz formula is used for more general permittivities as ([2]) 
is. Thereby it is always assumed that e(i£) is real such that the energy ((SJ) is 
real too. There are, however, also representations of Eq with integration over 
the real cu-axis, see Eq. (2.4) in [2] or Eq. (12.37) in [5], for example. These 
representations look like Eq. ([3]) with the w-integration starting from u = 
and no sums. It must be mentioned that these representations were derived 
under the assumption, motivated by dissipation or other physical reasons, that 
Im(e(u;)) > holds and that there are no poles of the transmission coefficients on 
the real cu-axis. Such systems do not have eigenvalues in terms of real frequencies 
oj in the sense as discussed above (see also the discussion in [B]). For this reason, 
there is no direct relation between such representations and Eq. ([3]). 

The Lifshitz formula is of fundamental importance in many areas of modern 
physics, including adhesion forces, atom-wall interactions, van der Waals forces 
and its applications reach to nano-technology and biology. Accordingly, this 
formula has a long history. Originally it was obtained in [2j from the energy 
of the electromagnetic field driven by fluctuating charges in the half-spaces, the 
latter being thermally averaged following the theory by Rytov. In [2J, also the 
limiting cases were derived. For large separation, Casimir's result was shown to 
follow. For small separation, a new formula for the van der Waals force, acting 
between half-spaces, emerged from (j5J) by putting in (jBJ) u p — > in x. The 
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corresponding interaction energy can be written in the form 
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(7) 




By rarefying the media in the half-spaces, also the known Casimir-Polder and 
the London forces between molecules emerged as limiting cases. 

It must be mentioned that the understanding of the Casimir and van der 
Waals forces appeared in [T] and in [2] from completely different ideas. While 
Casimir started from the zero-point energy of the electromagnetic field, Lifshitz 
considered the fluctuations in the medium as source. In the modern understand- 
ing these two are equivalent. However, the discussion about two ways continues 
until present time, see for example [7J. 

In the historical development, the next step 
after [2J was made in [8], by considering the 
problem in the general framework of thermal 
quantum field theory. This procedure was 
found very difficult and lengthy. Also, there 
was an interest in understanding the Lifshitz 
formula in terms of vacuum energy in the aim 
of Casimir. An important step in this di- 
rection was made by the very influential pa- 
per [9] in the non-retarded case starting from 
the vacuum energy ([T|). For small separation, 
when the electromagnetic interaction happens 
instantaneously, the dominating contribution 
to the mode sum in comes from the surface 
modes present in the TM polarization. From 
these, in j9], the limiting case ([7]) of the Lifshitz 
formula was obtained in a short derivation on 
2 pages. Soon after, this approach was gener- 
alized to the non-retarded case. In [TU], the idea in [H] of the vacuum energy 
as resulting from the surface modes was taken over literally to the non-retarded 
case by accounting for the full momentum dependence in the reflection coefficients 
OH]). The sum over these two modes in ([3]) was transformed into an integral over 
imaginary frequencies and formula fl5]) was reobtained. Independently, in [11] , 
this result was obtained by the same method, admitting, however, that, besides 
the surface modes, all evanescent modes should be included. 

It must be mentioned, that in the papers [10] and [11] the contributions from 
the photonic modes were missed. This was compensated by ignoring the cut 
in the complex frequency plane the functions generating the dispersion relations 
have. These omissions were observed quite easily, for example, in [T2] and in 




Figure 1: The configuration of 
two half-spaces with permittiv- 
ity s, separated by a gap of 
width L. 
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[13] (see also, the thesis [HJj, p. 76), but did not receive due attention. This is 
another reason for the present paper to clarify in detail the role of the various 
modes in ([3]). It must be mentioned, that the waveguide modes, including those 
in planar waveguides, are well known in electrodynamics, see [TH], for example. 
However, their role in the Casimir effect was, to the authors knowledge, never 
considered. One exception is the paper [16], where these modes were discussed 
(using different naming, see Eq. (8) there )q 

In the following section we consider the basic formulas for determining the 
modes and point out their basic properties. In Section 3 we consider the vacuum 
energy, separately for both polarizations. The last section has some conclusions. 
In the appendices we give a short description of the reverse transitions, from 
imaginary to real frequencies. 
In the following, we use units with h = c = 1. 

2 Basic equations and modes 

In this section we collect the basic formulas from classical electrodynamics and 
identify the modes. We do this in quite big detail. Although all these questions 
were considered and solved long ago, there is still some confusion about the ques- 
tion which modes are relevant for the Casimir effect in the given configuration. 

We consider two parallel half-spaces characterized by a permittivity e(uj) (see 
Fig. [I]). These are separated by an empty gap of width L with parallel interfaces 
and we have e — 1 in the gap. As it is well known, the solutions of the Maxwell 
equations for this configuration separate into TE and TM modes. After taking 
Fourier transform in the time t and in the directions X|| parallel to the interfaces, 
the electric field strength £ becomes proportional to 



where u is the frequency and k|| is the two-dimensional momentum in direction 
parallel to the interface. Since Eq. Q holds on the whole z-axis, we have to 
consider the permittivity e(u, z) with e(u>, z) = 1 for z G [0, L] and e(oj, z) = e(u) 
for z (fi [0, L], The components of the function <&(z) have to fulfill the well 
known matching conditions on the interfaces. For example, for the components 

1 This thesis can be ordered from the library of the University of Utrecht, the book- number 
is UB-ND MAG DISS UTRECHT QU 1975-25. 

2 I'm indebted to the authors of that paper for the hint. 




(8) 



and the Maxwell equations reduce to the equation 
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(10) 



corresponding to the electric field parallel to the interface these demand 
$(2) and £®( z ) (TEmode) 

and e(u,z)£$(z) (TM mode) 

to be continuous. 

The solutions of Eq. fl9]) can be written in the form 

{e ikz + R(k) e~ ikz {z < L) 

fie iqz + ue~ iqz (0<z<L) (11) 
T(k) e ikz (L < z) 

and represent a wave incoming from the left. Here, k has the meaning of the 
momentum in direction of the z-axis, i.e., in perpendicular to the interfaces, 
outside the gap, and q is the corresponding momentum inside the gap. In Eq. 
( fTTj) . R(k) and T(k) are the reflection and transmission coefficients, \x and v are 
some more coefficients. We indicate the momentum dependence of R(k) and T(k) 
keeping in mind that for the TM case there is an additional dependence on kn. 

From Eq. (Q, i.e., from the Maxwell equations, for these momenta the rela- 
tions 

e(co)uj 2 = k\ + k 2 , 

u 2 = k 2 + q 2 (12) 
follow. For the permittivity given by Eq. these relations read 



2 = ool + k^ + k 2 , 

2 _ !„2 , 2 



io 2 = k{ { + q 2 . (13) 

It is clear that out of the four quantities, u, k\\, k and q, only two are independent. 
Below we will use different choices of the independent ones, the remaining must 
be expressed in terms of these using the above equations. 

Using the matching conditions (flOi) . all coefficients in (TTT]) can be determined. 
In the following we need only the transmission coefficients, which take the well 
known form 



(q + k) 2 e-^ L - (q - k) 2 e^ L 



T TM (k) = Ae{u)qke * kL 

1 J (e(u)q + k) 2 e-^ - {e{u)q - k) 2 e^ L ' 1 ' 

for the two polarizations. For later use we rewrite these expression in the form 
T™(k) = ^W (£ ^f fc)a e^, (15) 
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with 



T 1 TE (/c) 



1 - r 2 E e i2qL 



where 
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k + q 



fc + e(u) q 

are the reflection coefficients on a single interface. 

By Eq. (jUJ) and the matching conditions ffTUj) . a spectral problem is set up. 
The spectrum consists of all real u, allowed by these equations and conditions. 
This problem can be illustrated by establishing a relation to a simple quantum 
mechanical Schrodinger equation by rewriting fl9]) in the form 

rf 2 „, A,,, 2 



dz 



2 +V(z))$(z)=u 2 $(z), (18) 



^ + fc 2 , (^[0,L]), 
= { 19 
^ (^[0,L]). 



with a potential 



With the matching conditions (flOl) for the TE mode, this is a simple exercise 
with a piecewise constant potential, the so called finite square well. For the TM 
mode the potential is more complicated. But anyway, we can determine the 
spectrum. In relation to the electromagnetic problem, given by Eq. (Q, it must 
be understood that we consider a problem with fixed but arbitrary k\\. 

For the TE mode, the solution of the spectral problem is obvious. We have 



scattering states with oj > y Up + and bound states with fc|| < u < y + 
the latter having imaginary momenta, 



k = iH = iJu^ + B-u 2 , (20) 



outside the box. As concerns the TM mode, it is known that it has the same 
properties (except for some bound states which can be located in < uj < ku 
too), although that is a bit more complicated to show. 

From the scattering problem, it is known that the transmission coefficients, 
T TE (/c) and T™(k), are meromorphic functions in the upper half of the complex 
fc-plane with simple poles on the upper part of the imaginary axis and a continu- 
ous continuation to the real A;-axis. The location of the poles is just the momenta 
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k = ix corresponding to bound states. For the function $(-2), this results in an 
exponential decrease for \z\ — > 00. Finally, we mention the property 



T(-k) = T*(k) 



(21) 



for real k under complex conjugation. It is this relation which motivates us to 
indicate the dependence on the argument k explicitly. 

Returning to electrodynamics, we can identify the following modes, 



The photonic modes correspond to the quantum mechanical scattering states 
and have a continuous spectrum. The waveguide and surface modes correspond 
to the bound states. Accounting for the momentum fen, these form a tower of 
continuous states. It must be mentioned that all electromagnetic waves in this 
problem are propagating waves having real frequency u. The difference between 
them is that the photonic modes do propagate in all spatial directions whereas 
the surface modes do propagate only in direction in parallel to the interface and 
the waveguide modes do propagate only inside the gap (including zig-zag like 
propagation). It should be mentioned that commonly waves with 00 < ku are 
called evanescent waves (see [5], p. 289). In Fig. [2] this corresponds to the region 
below the dotted line. 

Now we consider the waveguide modes in more detail. First, we take the TE 
case, which is easier since the transmission coefficient T TE (k) does not depend on 
k\\. As already mentioned, these modes correspond to the quantum mechanical 
bound states and their location is given by the poles of T TE (k). These appear 
for imaginary k = i-sjojp — q 2 - From Eqs. ({14]) . or (|T6|) together with (|T7|) . the 
conditions 



follow, which are for convenience written in terms of the variable q, see Eq. 
(|T3|) . The symmetry properties refer to the function $(2), (TTTT) . with respect to 
reflection on the middle of the gap. These conditions are exactly the same as 
that defining the mentioned bound states. The solutions of Eqs. ( 1221) will be 
denoted by qJ E with j = 1,2,..., [^] + 1, where [. . .] denotes the integer part. 
Using f lT3|) . we define the corresponding frequencies, 



• photonic modes with u > ,/io 2 + k?. 




• surface and waveguide modes with < u < ,/u| + E 




(22) 




(23) 



for which 




(24) 



S 



holds. The momenta k, related to these solutions, are imaginary, k — > ixj , with 



TE _ . / 2 _ („TE\2 



- (iD 2 - (25) 



The numbering is done in a way that 



uj E < w ™ (26) 

holds. Solving Eqs. (122]) numerically, pictures for the frequencies uJ E as function 
of fen can be generated. This is shown in Fig. [2J It must be mentioned that there 
is at least one such state for any fixed value of u p . 

A characteristic property of the waveguide modes is that these, for u p — >■ oo, 
turn into the modes known for ideally conducting walls, qJ E — > with j = 
1,2,.... This can be seen directly in Eq. (T22]) and this is equivalent to an infinite 
square well in place of ( |T9l) . 

Now we consider the TM modes. The condition for T™(k), Eq. ( fl4l) . to have 
poles, 

Jul - q 2 { tan^ (symmetric) 

V = < (27) 

£ { U} )Q ^ — cot ^ (antisymmetric), 

differs from (1221) by the presence of the permittivity, 

q 2 + kl - u)l 

only, where (|T3j) was used to express a; in terms of q. We denote the solutions of 
(|27|) by gJ M . For fcy > w p we see immediately that, because of e(u>) < 1 in the 
left hand side of (12 7p . the TM solutions are larger than the corresponding TE 
solutions, 

qJ M >qJ E , (J = l,2,...,pl+1). (29) 

L 7T J 

By means of ( 1131) . a similar relation holds for the frequencies, 

^™ > coJ E , (30) 



where we defined 



oo 



For fc|| — >■ oo, we have e(o;) — > 1 and both kinds of solutions coincide. In this way, 
we observe a one-to-one correspondence between the TE and the TM modes. For 
fen < u p , the picture changes since e(u) as a function of q goes through zero. In 
this case the mentioned correspondence persist, but an additional solution, q™ 
with a corresponding frequency u™ (also obeying the inequality (l2"6"j) ). appears 



provided hi < uJ y 1 + u p L/2 holds. 
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Figure 2: The spectrum of the electromagnetic field for uj p = l.lix. Solid lines rep- 
resent, from bottom to top, the frequencies of the symmetric and of the antisym- 
metric surface modes and, further, of the TM waveguide modes uj™ (j = 1, 2, 3). 
The dashed lines are the frequencies of the TE waveguide modes ujJ e (j — 1, 2, 3). 
The dotted line corresponds toco = k\\ and separates imaginary {oj < k\\) from real 

values of the momentum q (see Eq. ( fl2|) ). The shaded region has uj > 

and represents the photonic modes. 



Now, since all q™ are real, for the frequencies the inequalities 



fc„ < cJ M < Jufi p + B 



(32) 



TM 



hold like Eq. f )24l) in the TE case. Also for these modes, the momenta k — > ixj 
with 

(33) 



TM 



uj- 



TM\2 



are imaginary. 

There is one more relation between the TE and the TM modes. For k\\ — > 0, 
Eq. (|2T|) for the TM modes turns into Eq. (|22|) for the TE modes. Thereby the 
correspondence changes by one. In this way, the inequality 



, ,TM ^ , TE < TM 



(J = 1,2, 



UJ,, 



71 



(34) 
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can be established. The left equality is reached for k\\ — y and the right one for 
fc|| —> oo. 

In addition to the TM modes considered above, there are also poles of the 
reflection coefficient for imaginary momentum, q — y ir], for which Eqs. (1271) can 
be rewritten in the form 



y/u* + rj 1 ( tanh^ (symmetric), 
e(u)r] y coth ^ (antisymmetric), 
and the permittivity expressed in terms of i] reads 



(35) 



£ ^ = k 2_ v 2 ■ ( 36 ) 

Each of these equations has one solution which we denote by r) s f (sf = s,a 
according to symmetry). The corresponding frequencies are 



u*f = \J k l- V 2 sf ■ (37) 
Of course, also in this case the momenta k are imaginary, k — y ix, with 

"•f = y/ u P + tfr ( 38 ) 

These two solutions, co s and u a , correspond to waves decreasing exponentially to 
each side of each interface in opposite to the waveguide modes which decrease 
only outside the gap but oscillate inside. 

It is known that, for L — y oo, these solutions turn into the surface plasmons 
well known in the physics of metals. In this limit, the right hand side of Eq. (j35p 
equals unity and the equation allows for an explicit solution, 



^single 



4 + K-\(4)' + K- (39) 



At large separation, such solution exists on each interface. At finite separation, 
the frequency 0J sing i e splits into two, co s and co a , for which the inequality 

W s < ^single < (40) 

holds. 

There is one peculiarity about the antisymmetric surface plasmon. It is a 
solution of Eq. f )35|) only if the inequality 

< h (41) 
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holds. In that case the momentum in the gap, q = ir] a} is imaginary. For smaller 
feii, this mode matches the solution u™ which was found among the waveguide 
modes if just the opposite to (I4TT) holds. It is clear that these constitute one single 
mode which must be identified with the antisymmetric surface mode. We mention 
that the question of whether to include the lower part of the antisymmetric 
surface plasmon, where its frequency is given by u™, into the vacuum energy, 
was discussed in [T7] and [IS] . 

The frequencies of the surface modes are shown in Fig. [21 too. The antisym- 
metric one crosses the line u = fen just for fen = tu p / a/1 + co p L/2. For lower k\\, is 
is given by co™, Eq. ( IBTj) . and, for larger fen, by u a , Eq. ([3"Tj) . 

Finally we mention, as a quite special property of the surface modes, that 
their frequencies have an upper bound, 

Usf ~ % ■ (42) 

This can be seen by mentioning that u s f are monoton functions increasing with 
increasing fe||. For large fe||, the solutions 77 of the equations (1531) grow proportional 
to fen. Consequently, the hyperbolic functions in the right hand side of these 
equations turn into unity and we are left with the equation for a single surface 
plasmon. From the corresponding energy, Eq. ( 1391) . one can then infer the bound 



3 Vacuum energy as mode sum 

In this section we consider the vacuum energy of the electromagnetic field in 
the sense of Eq. ([1]) and will give Eq. fl3]) a precise meaning. We start with 
mentioning that the spectrum, if accounting for k\\, is completely continuous. 
Hence one needs to separate the empty space contribution. This can be done 
in numerous different approaches. We follow that used in [T9]. There, a large 
box was introduced, whose volume was subsequently tended to infinity. As a 
result, Eq. ([3]) appears with the scattering phase shift expressed in terms of the 
transmission coefficient, 

S = — In — ^-f . 43 

2i T(k)* v ; 

Another, equivalent, form of writing is 

8 = Im(lnT(fe)), (44) 

which is also frequently used in literature. 

Because of the separation of polarizations in the considered problem, the 
energy consist of two parts, 

E = E TE + E™. (45) 
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In the next subsection we consider the contribution from the TE polarization. It 
is easier to handle. In a subsequent subsection, we consider the TM case. 



3.1 TE case 

For the TE polarization we have to consider waveguide and photonic modes. 
Accordingly, we divide the regularized vacuum energy, 

E TE = + Ej E t} (46) 

into the waveguide contribution, 

^H/^EW + l^) 2 )*"*. (47) 



J 



and the contribution from the photonic modes, 

1 r dkn f°° dk . 2 , 2 , 2N i- s d , T TE (k) , , 

As compared to Eq.Q, we changed the variable of integration from u to k us- 
ing (fT3|) . It must be mentioned that both above expressions carry ultraviolet 
divergencies and one needs to keep the regularization with s > |. 

Now we carry out the integration over k||, which can be done in the TE case 
since neither T TE (k) nor qJ E depend on k|| (note q = + k 2 from Eq. (IT3"])). 
The integration is simple and we come to 

E TE = — V (a TE ) 3 ~ 2s (49) 



and 



^ = s^r^M + *') , -s h ^©- (50) 

In the waveguide contribution -E™; ^ ne analytic continuation can now be carried 
out simply by putting s — 0. However, we postpone this. The second contribution 
will be split, using the factorization in (fl5j) . 

= £ Jont + #£l> ( 51 ) 



into 



and 



E ^^^ ) [^^ + ^ i "^- k) - (53) 
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Note that the factor 4qk/ (k + q) 2 dropped out since k and q are real here. 

The latter contribution is proportional to the width L of the gap and it can 
be calculated explicitly, 

2(2-s) j 

E Ll = -^^h 1 {sl (54) 

with 

1 v/?r(s-2 
2(s^2) ~ 2r(s - f ) 

The function h\(s) has a pole in s = 0, 



Ms) = 7771 n\ - ^T^V- (55) 



h 1 (s) = -^ + 0(l), (56) 
lbs 

which results from ultraviolet divergence. 

In the first contribution, E™ t , it is possible to put s = 0, 

^a=^jff M+^'l^w), (57) 

since the integral is convergent due to the decrease, 

l nTl T Em^4ie m for Awoc, (58) 
16k 4 

which can be checked using ( 1T5|) . 

Taken in the form of Eq. fl52l) . we can calculate -E™ t numerically since it is 
represented by a convergent integral. For technical purposes, it is useful to turn 
the integration path slightly up into the complex plane, k — > ke ia with a > 0. 
Under such substitution -E™t does not change, but the integral converges much 
more rapidly. 

Next we establish the relation to the representation on the imaginary fre- 
quency axis in order to get the TE part of the Lifshitz formula, (j3j). We start 
from (1501) . which we rewrite in the form 

E TE - - 1 r dk U + k 2 ) ^ d In TTE{k) (59) 
^ cont " 4vr(3 - 2a) J Q 2m K + k > dk^T^k)' (59) 

using (|2T|) and keep s > |. We are going to change the integration path in ( 1501) 
following the procedure in [19] . The complex fc-plane is shown in Fig. [3j The 
integration in fl59|) is over the real A;- axis, k e [0, oo). On the imaginary axis, the 
integrand has simple poles with residua equal to (-1) at the locations k = ixJ E , 
Eq. ( }2"5|) . Above, starting from k = iu p , a cut starts which results from the factor 

(uj 2 + A; 2 ) 3 ^ 2 s in the integrand. Furthermore we mention that T TE (k) has a zero 
in k = 0. In the quotient in the logarithm in ( l5"9"j) these zeros cancel and the 
integrand is regular in k = 0. 
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Figure 3: The complex &;-plane. The 
dots represent the poles of the in- 
tegrands in (jnTJ and (I106p resulting 
from zeros of the reflection coefficients, 
the crosses represent poles resulting 
from the poles of the reflection coeffi- 
cients. Their locations correspond to 
the waveguide and surface modes whose 
frequencies are shown in Fig. |2j The 
pathes 7 r and 7 are used in Eqs. (1ST]) 
and (I106p . On the imaginary axis there 
is a cut, starting from u p in the TE case 

(see Eq. (1651) or from 
case (see Eq. (1109]) ). 



u 2 p + k\ (TM 



Im(k) 



Re(k) 



Now we represent the logarithm in floTT) as a difference, lnT TE (fc)— lnT TE (— k), 
and split the integral into two accordingly. In the second integral we change the 
variable, k — )■ —k. In doing this we pay attention to the pole in k = 0. Therefore 
we write 



E. 



TE 
cont 



1 



, lim 

4tt(3 - 2s) r^O 



+ 



dk 
2hr 



d 



K + fc 2 r S ^lnT TE (fc), (60) 



It is our intention to move the integration path upward in the complex plane. For 
this, we need to close the path by adding and subtracting a path half encircling 
the pole in k = 0, 



E 



TE 

cont 



-1 



, lim 

4tt(3 - 2s) r->o 



* K + ^.„T- W . (61) 



The paths 7 T and 7 are shown in Fig. [3j In the limit r — > 0, the contribution 
from the path 7 T picks up half a pole in k = and the path 7 becomes independent 
of r. We get 



pTE 
-^cont 



-1 



4vr(3 - 2s) 



1 , ?- 2s 
2 W p 



dk 



(62) 



Now, we move the path 7 upward to k 



We cross the poles in k = Ik, 



TE 



which delivers just the sum S™, Eq. ( 1491) . with opposite sign. After that we 
tighten the path around the cut on the imaginary axis, k — > ix, and with 



( Wp 2 + M 2 )- s -h 2 + H,f) 



-2i cos(7rs) (x 2 



4Y~ 



we get 



pTE 
-^cont 



4tt(3 - 2s) 



x , ,3-2s 771TE 



+ E 



TE 



(63) 
(64) 



15 



where 



is the unrenormalized vacuum energy represented as integral over the imaginary 
axis. It is, however, still different from the corresponding expression in ([5]) in 
having the transmission coefficient T TE instead of T™ in ([5]). To proceed, we 
make use of the factorization ( TT51) which, analytically continued to the imaginary 
axis, takes the form 

T TE (*x) = TfV) , 4r/X , 9 e-b-** (66) 



For the analytic continuation of the momenta we use the notations k = ix, q = %r\. 
The relation q 2 = uj 2 + k 2 , following from (JTHjl . translates into r\ 
According to (1561) we split 



2 - x 2 - 



^imlg = ^hnag + #TE + E Li , (67) 



where, from the first factor in 



from the second factor 



cos(7rs) Z" 00 , / o o\ S-s <^ , 4r?x , , 

f ™ = mH5JL*' ( ^" 6 * (69) 

and from the third factor 

_, C0S(7TS)L f 00 , / 2 2\%~a & , N 

^ = -4^4)4^^-^ ^ (T <- X) (70) 

follow. The latter contribution is proportional to the separation L and can be 
calculated easily, 

2(2-s) r 

^ = - ! W £ W')- < 71 > 

The function ^(s) can be expressed in terms of hypergeometric functions. It has 
a pole in s = 0, 

h 2 (s) = ^ + 0(l), (72) 
16s 

with the same residuum as hi(s), Eq. fl56l) . In the following, we will need only 
the relation 

\i m (h 1 (s)-h 2 (s)) = -\. (73) 
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Now we consider , Eq. ( )68|) . First, we note that the function 

with r TE given by Eq. (E}, is just the same as in the first logarithm in (j5j. 
Since it is exponentially decreasing for x — > oo, we can remove the regularization 
by simply putting s = 0. Finally, we perform a substitution of variables, £ = 

\/ x2 - w p > 

T™{i*) = -\ 2 (76) 

1 - {U ) e " 2eL 

is the transmission coefficient on the imaginary axis appearing in (J5J) after carrying 
out the integration over k\\ . Expression ( 1751) can be considered final. Because of 
the convergence, it is also possible to integrate by parts, 



where 



E, 



TE 

imag 4?r 2 



l - 2 [ 'dteinT™(ix). (77) 



Eqs. (1731) and (177j) are variants of representing the well-known TE part of the 
vacuum energy calculated on the imaginary frequency axis. We mention that 
these differ from E^ ag , which we initially got from moving the integration path, 
by the addendum H TE , which does not depend on the separation L, and by 
El 2 , which is proportional to L. In the sense of a renormalization it is natural to 
remove these two contributions and to consider as the renormalized vacuum 
energy. Of course, -E™ g has the necessary behavior for large separation, 

E ™^-xm? for L ^°°- (78) 

It is also easy to obtain the behavior for small separation. One can put L = in 
the integrand in (1751) since the remaining integral converges. A short calculation 
gives 

3tt - 8 

727T 2 

This limit coincides with that found in [19], Eq. (41). In the complete energy, 
this contribution is subleading. The leading contributions comes from the TM 
polarization, Eq. (j!2ip . 

It remains to consider i^TE ; Eq. ( 169]) . The integration can be carried out 
explicitly, 

COs(7Ts) f 1 4 



E Z e ^-^-r^ r for L-,0. (79) 



"-i^lrsl"" < 80 > 
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With this, we return to E™ t , Eq. (16"4"|) . and, inserting from (157)1 and ( 180|) . we get 

7~>XE COS^7Tt5^ 3— 2s r~iTE rn /oi N 

-^cont = ~3(3 _ 2s)7T 2 Wp + ima § + L2 ' ' ' 

The contribution from the half pole in = canceled. We mention that we have 
to keep s > | (in fact, s > is sufficient) in this expression because of the pole 
in E L2 . 

The equations (IHll and ()46l) . together with (1471) and fl5TT) . are two represen- 
tations for E TE , one with integration on the imaginary axis, and the other with 
summation and integration on the real axis, 



771TE _ V-/ 1 piB 771 

E ~ "3(3-2*)^ + £ w g + ^ 2 

= -^wg + + ^Xj • ( 82 ) 

Needless to stress that these are equal, at least for s > §. Now we perform 
the renormalization by subtracting from both sides the first two contributions 
in the first line. The first of these does not depend on L and the other, El 2 , is 
proportional to L. Thus we define 

77TE _ g.TE 1 COs(7Tg) 3 _ 2s 

E ~ E + 3(3- 2 a )7r> W * ~ El > m 

as the renormalized vacuum energy. It has all necessary properties since from 
fl82D and flU 

E TE = E™ g (84) 

follows by definition, for which these properties are known. 

Now we insert the second line from (1521) into (1531) and come to 



COS! TVS I 



E TE = £v g + £™ t , (85) 
where -E™ t is given by Eq. (1521) . and where we defined 

P TE _ pTE 1 CQS(7TS) 3 _ 2s 



Now we are finally in a position to remove the regularization. The pole in s = 
is present only in the last two terms in ( 15B"|) and cancels. Using (173"|) together 
with ( )54|) and (1711). we get for s = 



pTE 1 \p /TE\3 U l /o 7 N 



12tt v ^ y 9vr 2 48vr 2 



3 



We included the contribution proportional to L, resulting from E Ll and E L2 , and 
the constant contribution, resulting from H, into the waveguide contribution for 
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Figure 4: The contributions to the vacuum energy for the TE polarization as 
function of u p for L — 1. The solid line represents the vacuum energy, -E™ g = 
+ E™ t , the thin solid line represents the contribution from the photonic 
modes, E™ t , and the dashed line is the contribution from the waveguide modes, 
E^. The curves are equivalent to represent the energies multiplied by L 3 as 
function of L for u p — 1. The jumps appear when u p goes through multiples of 

71. 



the following reason. It can be shown, that the sum in (1BTI) has for large L an 
asymptotic behavior 

E ™* = ~48^~9^ + ---' (88) 
where the dots denote bounded, but oscillating contributions. 

In (IHTj) . the first two terms are just subtracted. Hence, E™, Eq. ([57]) . at 
most oscillates for L — > oo, which justifies considering it as the contribution from 
the waveguide modes to the vacuum energy after renormalization. 

Collecting from f)84p and 0851) . we have with 

771 TE _ T-i TE _ ttiTE j_ rnTE /o Q \ 

two representations of the renormalized vacuum energy resulting from the TE 
polarization, one with integration over the imaginary axis, which coincides with 
the corresponding part in the Lifshitz formula, and another one with summation 
and integration over real frequencies, i.e., over the spectrum. Both sides can be 
evaluated numerically. Both must deliver the same numbers. Indeed, they do. 
The results are represented in Fig. H]as a function of u p with L — 1. It must be 
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mentioned that, already for dimensional reasons, the energy can be represented 
in the form 

E TE = ±f(u p L), (90) 

where f(u p L) is dimensionless (its limiting values can be read off from Eqs. (j78p 
and fl79l) ). Thus a plot of E TE as function of u p is equivalent to a plot of E TE L 3 
as function of L. 

In Fig. H] also E^ E and -E™t are shown as functions of u p . These show an 
oscillating behavior with strong compensation in a way that their sum is just the 
much smaller , shown as thick solid line. 



3.2 TM case 

In this subsection we consider the TM part of the vacuum energy. We start from 
the unrenormalized one, E , in Eq. (1451) . According to the spectrum of modes 
described in Section 2, it is a sum, 

E™ = E si + E™ + E™, (91) 



of surface modes, 



of waveguide modes, 



wg ~2J (27p Z^\ u i > ' (93 ^ 
and of photonic modes, 

~ TM _ If dk„ Sk f°° dk ( 2 + k 2 + k2) hs d T™(fc) 

It must be mentioned that we are forced to introduce an additional regularization 
for the integration over the momentum k||. The reason is in the boundedness of 
the frequencies of the surface modes, as expressed by Eq. fT4"2l . which makes the 
zetafunctional regularization effectless. This is some kind of additional divergence 
which was not observed earlier. However, it does not show up in the force. As 
we will see below, all terms divergent for 5 — > 0, do not depend on the width L 
of the gap. 

The main problem, we are faced with, is to perform the analytic continuation 
in the zetafunctional parameter s to s = 0. In the TE case this was simple 
since, after performing the integration over k||, we could perform the analytic 
continuation explicitly. In the TM case we cannot integrate over kii in such 
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simple way since it enters the reflection coefficient. However, below we will show 
how this problem can be solved. 

We start with the surface mode contribution ( J92l) . Here we can make use 
of the circumstance discussed at the end of Section 2, that the frequencies, for 
large fcii, approach the frequency oj single-, ( )39l) . of a surface plasmon on a single 
interface. Therefore we subtract u single and add it back, 



where 



E s i = E si + 2E sing i e , (95) 



sf 

is the subtracted energy of the surface modes and 

E single = ~ J ^Jy *S <5fc| 1 ^ (^single) 1 ^ (97) 

is the energy of a plasmon on a single interface. In the energy E s f one can remove 
the regularizations, i.e., one can put s = and 5 = 0, and the emerging integral 
turns out to be convergent. The integral in E sing i e , Eq. ( 197|) is quite simple and 
can be calculated explicitly for s = and 5 — > 0, 

£„„„. = ( 1 " 12 '" 2 ± 6 '""' + M^±2l ) 4 + la + (f), (98) 
V 192^ 32^ / p 4v^7r<J 2 

where 7 is the Euler constant. It has two terms divergent for S — > 0, one quadratic 
and one logarithmic. We mention that E S i ng i e does not depend (by construction) 
on the separation L and therefore all these terms depend on u p only. 

Next we turn to the waveguide mode contribution, Eq. (1931) . First, we mention 
that we can remove there the additional regularization putting 5 = 0. This is 
because the frequencies u™ of the waveguide modes grow proportional to k\\ 
for fc|| — > 00 and the zetafunctional regularization is effective here. In order to 
perform the analytic continuation in s, we subtract the frequencies of the TE 
waveguide modes and add them back, 



E™ = AE wg + E 1 ™, (99) 



where 



has the differences of the frequencies and i?™ is the TE contribution (|47p . The 
latter has been treated in the preceding subsection and the analytic continuation 
in s is given by Eq, (j4"9~|) . The difference energy, AE wg , can be continued in s 
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simply by putting s = in fllOOp since both frequencies approach each other for 
fey — » oo, as mentioned after Eq. (l3Tj) . 

Finally, we consider the energy E™, flE}, of the photonic modes in (I9TT) . We 
proceed similar as just above and subtract the corresponding TE frequencies, 

= A-^cont + £™ t , (101) 

where E™ t is given by Eq. p5|l and its analytic continuation by the subsequent 
equations in the preceding subsection. The difference energy, 



d 

dk 



(T™(fc))* (T TE (fc))*. 



(102) 



has the property that on can put s = directly since the remaining integrals do 
converge. Also we can put 5 = since dangerous terms cancel in the ratios inside 
the logarithms. For numerical evaluation we used 

i noo roo id/ T™(k)\ 

AE cont = _ j( dk ^ d*„*|| (ul + % + Im [\n _Uj , (103) 

which differs from f ll02p by simple rewriting, including a change in the order of 
the integrations. 

Now, after having clarified how to make the analytic continuation in the 
contributions from the real frequencies, we turn to the imaginary frequencies. 
Again, we follow the procedure used in the preceding subsection. We are going 
to transform the integration path of the k integration. We do this for fixed k\\. 
For that reason, it is useful to start from the contribution of the photonic modes, 
E™ t , (194p . which we represent in the form 

EZ=lJ^e-^F(k n ), (104) 

with, using (12~TT) . 

r°° dk id T™(k} 

m)=[ £K+^«T-s to TO- (105) 

The transmission coefficient T™(k) has properties similar to T TE (fc), except for 
an additional zero for k = iku. For instance, it has a zero in k — 0. So, from 
the ratio inside the logarithm, we split the logarithm as a difference of two and, 
further, the integral as a difference of two with substitution k — > — k in the second. 
Accounting for the pole of the integrand in k = we come to 
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We proceed by closing the integration path around k = as explained in the TE 
case resulting in 

F W=Ss(/ + /J^M + *') , "*s h ' I ™w- (107) 

Now, as mentioned, we have an additional pole in k = iku. It may be located 
inside the contour 7 r or outside in dependence on k\\. It can be seen, that the 
case, when this pole is located inside, results in a contribution to (1104p having an 
integration over k|| with k\\ < r which has for r — > an additional smallness ~ r 2 
and the corresponding contribution vanishes in the limit. So we are left with the 
case that this pole is outside. Then the contribution from 7 r can be calculated 
directly and, as before, 7 becomes independent on r, and we get 



nh) = \ H + *5) + H + *S + * ! ) ,_ * Tk lnT ™ w ■ (108) 



Now, we move the path upward until k = i\/Up + k?,, crossing poles of the inte- 



grand. These result from the poles of T (k) in the points k = ix s f (correspond- 
ing to the surface modes) and from the poles of T (k) in the points k = %k. 



.1 



(corresponding to the waveguide modes) having residua with negative sign. Fur- 
ther we have a pole from the zero of T™(k) in k = ik\\ (which was not present 
in the TE case). Finally, we tighten the path around the cut starting now from 

k = i./Up + k?, and come to the representation 



sf j 



+ cos(7T£) r° rfx ( x 2_ uj.fcjj)*- A. lnT™(zx). (109) 

We mention that the cosine appears with opposite sign as compared to (|63|) 
because of the different number in front of s in the exponential. 

Now we insert (I109p into (11041) . The integration over k\\ can be carried out in 
the first two contributions, for the next two we use fl92l) and ([93]) and come to 



pTM _ U P , 1 U P B pTM , pTM / lln \ 

C ° nt ~ _ 247 4^2 52 ~~ hs f ~ ^wg + Amag, I 11 ") 



where we defined 



~ TM _ COs(Trs) f dk\\ Sk f°° , 2 , 2) l-s d , 



TM/ • \ 



111) 
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The latter is the unrenormalized vacuum energy represented as an integral over 
the imaginary axis. It differs from the corresponding contribution in (JSJ) by the 
different transmission coefficients, T (ix) here and T™(ix), Eq. (fT6l) . in ([5]). 

Now we consider this difference using (1T51) . analytically continued to the imag- 
inary axis, where it takes the form 

t tm m = t tm ( ,^ ) JjjtgBLp e~^. (112) 
We use the same notations as in Eq. ( 16T)j) and, in addition, the analytic continu- 



ation, to = it;, of the frequency, resulting in the relations £ = yx 2 — — /cjj and 

s(i£) = 1 + |f = Jf2 _ fc 2_ faJ 2 > which must be used in ( 11121) . 
According to the factors in (I112p . we split, 

E ™ g - E™ s + Htm + E L2 , (113) 

where El 2 is the same as in the TE case, Eq. ( 1701) with its analytic continua- 
tion ( TTTT) . In the first contribution, we can remove the regularization since the 
integrations are convergent such that it takes the form 

E ™ 8 =i\wff^"^- * ~ k ^i ln r ™ M - (U4) 

The reflection coefficient on the imaginary axis is 

TTM = ; i—^z, (H5) 

1 - r TM e V 

with ttm given by Eq, (J6]). Thus, Eq. (j!14p is just the TM contribution in 
the Lifshitz formula (151). It can be rewritten also by making a change of the 
integration variable x for £, 

or, by integrating by parts, in the form 



£™ = ^ / </M- / ^ ln:/?' M i^). (117) 



oo 



In (11161) and (I117p the variables must be expressed according to 



x = Jf + ut + M, 



-2 



e^O = 1 + ^- (118) 
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For later use we define the difference between the corresponding TM and TE 
contributions, 

A rp _ pTM _ rpTE (-\-\Q'\ 
'-'"imag — "imag "imag- l J -- Li V 

The TM contributions has the following asymptotic properties. For large separa- 
tion L, or, equivalently, for large u p , it turns into the ideal conductor expression, 

£™ g ^-14i5Z3 for ^">°°> ( 12 °) 

in parallel to fl75|) . For small separation, or small u p , we make in (11171) the 
substitution £ — > After that, one can put u p = in the integrand and come 
to 

£™ g ( 121 ) 

where the remaining integral, 

1 r°° r°° ( e - 2k w V 1 

C ^*W» *"H ^{'-(TTWr) • (122) 

is a number, C2 — 0.00391, which coincides with that found in [19], Eq. (41) 
and earlier in Eq. (14) in [20] and in [21 j . where a sum representation was 
derived. This expression can also be calculated starting from the contribution 
of the surface modes, Eq. (I147p . The corresponding calculation is shown in 
Appendix B. We mention that Eq. f lmj) with ffT22|) is the same as with 
inserted. 

Finally, the contribution 

in (I113p . resulting from the middle factor in (I112p . needs to be considered. In 
(11231) we changed the integration over x for £ ( fll!8j) must be used) and introduced 
the notation 

#11,0 = 1.^. (124) 
To proceed, we represent this function as a sum of three, 

= 9a + 9b + 9te, (125) 

where 

9A ~ ln w 

4(o; p 2 + 2e 2 )(x + r / ) 
or = 2 m ; — — — , 

9te = In t — ■ — r=. (126) 
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According to the splitting (11251) . we represent (11231) as a sum, 

H TM = H A + H B + H TE (127) 

The last term, following from the function (?te> is the same as appeared in the TE 
case, Eq. (1691) and can be treated in the same way as in the preceding subsection. 
The result is given by Eq. (|80l . 

Now we consider Ha- Since g A depends on £ only, the integrations can be 
carried out easily. For s = and 5 — > we get 

H ^^ 1 ^ + 0(5) - (128) 

This term has a quadratic divergence in the regularization parameter 6. This is, 
after E s f (see Eqs. (|95|) and (|98|) ). the second place where it appears. 
Finally, we consider H B , 

with given by the second line in (11261) . In order to perform the continuation to 
s = and to 5 = 0, we subtract and add back the asymptotics of g B for ku — > oo, 

9B ~ 25^+2£%|+if- (130) 
splitting into two parts, 

H B = H Bl +H B2 . (131) 
In the first part, -Hb i; the integrations converge for s = and 5 = 0, 

J /» (J\<i poo J I | 

H ^ = ^jwfL ^rww 1, (132) 



These integrations can be carried out numerically, resulting in 

c 
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H b, = — 2 "1 (133) 



with c ~ 0.06777. In the second part, 

i r dk\\ r°° d l 

HB2 = 2^Jj2^J da ^2(^ + 2e)K + A;2)' (13l) 
the integration can be carried out analytically in the limit 5 — > for s = 0, 

HB * = {H ieX'' ) "' + 0{S) ' (135) 

26 



where 7 is the Euler constant. Together, from f |133[) and fl 1 3 5 j) . we get 

«e=(^ + ^£^U + {S) . (136) 



.2tt 2 16^ 
This must be inserted, together with d 128[) and (1501) . into (11271) . 

^™-^-9^ + 2^ + l6V2n ) UJp + ^r¥ + ° iSh (13?) 

which completes the analytic continuation in E™ g , Eq. ( 11130 . Now we insert 
this expression into E™ t , Eq. (11101) and obtain 

pTM _ p pTM, ( C 1 1 (jgHj)+7) V,3 I pTM , 77 1 ,1 ooN 
^cont - ~^sf - ^ wg + - ^2 + J W P + ^mag + E L 2 , l 1 ^ 8 ) 

whereby four terms canceled. The last step in this sequence of calculations is to 
insert this expression for E™ t into (IPTT) . 

E -^"9^ + ^oV!^J p+ imas+ L2 ' (139) 

whereby the contributions from the surface and waveguide modes canceled. 

With Eq. (I139p . the transition from the summations and integration over real 
frequencies to the integration over imaginary frequencies, including the analytic 
continuations, is finished. Now we can define the renormalized vacuum energy 
resulting from the TM polarization. We proceed similar to the TE case and define 

-.TM _ 77TM / C 

1 (ln(w p 5)+7)\ 3 



EL ™= ElM -{^-^ + "-^£M^r'- E ^ (140) 

This coincides, by definition, with the vacuum energy -E™ g on the imaginary 
axis, 

E™ = (141) 

As already mentioned, E™ g is just the TM contribution to the Lifshitz formula 
©• 

In order to get the corresponding expression on the real axis, we insert for 
E™ in ( 11401) from Qj2). We get 

E™ = E™ + E™ + E™, (142) 

where we defined 

-^cont = -^cont ~ E L X - (143) 

Using (11011) and (15T1) . this can be rewritten as, 

E™=Ej l t + A£ cont , (144) 
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and represents the renormalized contribution from the TM photonic modes. It 
is represented as sum of the corresponding TE modes and the difference between 
TM and TE mode's contribution. 

The first contribution in E , Eq. ( 11421) . is defined as 



Kf = Ki ~ ~ { + ) ^ ' (145) 



and the second is defined as 



-.tm 771TM ■ ' c ' 1 12 In 2 



The terms, which are subtracted in these two expressions, follow from the defi- 
nition (I140p of the renormalized vacuum energy. Thereby the definition ( I145j) is 
motivated from Eq. (|98p in the sense, that just the doubled contribution from a 
plasmon on a single surface is subtracted and 

E S J = E s f — 2E ' S i n gl e (147) 

holds. After that, the subtractions left for (I146p are already uniquely determined. 
These can be shown to be just that substraction, which makes E™ vanishing for 
large separation, L — > 00, thus making E™ the correctly renormalized contri- 
bution to the vacuum energy. The waveguide contribution can be also split into 
the TE contribution and a remainder, 

E™ = E™ + AE wg , (148) 

where E^ g was defined in (157)1 and AE wg is 

( c 1-I21n2\ , . , 

with AE wg defined in fllOOf) taken with s = and (J73l) was used. 

Collecting from the above formulas, we have a with (11411) and (11421) two 
representations of the vacuums energy resulting from the TM polarization, 

7-.TM _ 7-.TM _ tttTM 1 771TM 1 P TM /-. r n \ 

E - Amag - E sf + -& w g +-& cont . (15U) 

One is with integration over imaginary frequencies, the other is over real fre- 
quencies, i.e., over the real spectrum. Like the TE case, both can be evaluated 
numerically and must deliver the same numbers. Indeed, we checked, they do. 
For technical reasons we represent the results as difference between the TM and 
TE cases 

AE = E™ - E TE . (151) 
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Figure 5: The difference between the contributions to the vacuum energies re- 
sulting from TM and TE polarizations from the photonic modes (solid line) and 
from the waveguide modes (dashed line) as function of ui p for L = 1. The dotted 
line represents the sum of the energies of the surface plasmons. The thick solid 
line is the difference between the complete contributions from the TM and the 
TE polarizations. 



Using ( TTT9D . (fl44l) and (TLTHl) . together with (EH]), we get 

AE = A£ imag = E™ + AE wg + AE cont (152) 

The surface modes appear in the TM part only, therefore these remain without 
TE counterpart which could be subtracted. The constituent energies in Eq. ( 11 52f) 
are shown in Fig. |5]as functions of the plasma frequency u p . Like in the TE case 
in Fig. HJ there are oscillations and compensation. All curves decrease for large 
uj p showing that the difference between TE and TM becomes subleading for large 
separation or large plasma frequency. 



4 Conclusions 

In this paper we have discussed the role of the electromagnetic modes contribu- 
tions to the vacuum energy for two half-spaces characterized by a permittivity e 
and separated by a gap of width L with e = 1. In this configuration, the renor- 
malized vacuum energy can be represented as sum of contributions from the TE 
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and the TM polarizations, 



E = E + E , (153) 

where E TE is given either in terms of integration over imaginary frequency, ( 1841) 
for E TE and (11171) for E™, or in terms of real frequencies, ( 1891) for E TE and (11421) 
for E™. Of course, both kinds of representations are equivalent, and for both 
we derived representations suitable for numerical evaluation. The corresponding 
numbers have been checked to coincide within the confines of reasonable numer- 
ical precision. Results are shown in Figures H] and [5j 

The representation in terms of imaginary frequency is, of course, identical to 
the Lifshitz formula, Eq. ([5]). As discussed in the Introduction, the representa- 
tion in terms of real frequencies, considered here for the plasma model, cannot 
be compared directly with those representations assuming Im e(cu) > which 
can be found in literature. Restricted to the plasma model considered in this 
paper, it is important to stress that the spectrum, i.e., the set of contributing 
real frequencies, consists of surface (for TM polarization), waveguide and pho- 
tonic modes. While the surface modes are well known to contribute, the role of 
waveguide and photonic modes is less clear in literature. As already mentioned 
in the Introduction, in pU] only the surface modes were considered, in [IT], in 
addition, all evanescent waves were included. Also in [I], it was not observed 
that there is a difference between waveguide and photonic modes. However, in 
[T6] attention was paid to this difference. However, that remained without conse- 
quences for the calculations and conclusions of the mentioned papers (and many 
others) since all calculations were performed in fact on the imaginary frequency 
axis and only the contribution from the surface plasmons were calculated on the 
real frequency axis. The photonic and the waveguide contributions were never 
calculated separately for the plasma modelE 

In Section 3, with Eqs. (189]) and (11421) . we derived representations of the 
vacuum energy in terms of real frequencies. The renormalization was done in a 
way as to ensure decrease for L — >■ oo. This was archived, in the TE part, by Eq. 
(I53"j) . and, for the TM part, by Eq. (I140p . In each case, the subtraction has a term 
which is independent on the separation and one term which is proportional to 
L. The latter appears as part of the subtraction of the empty space contribution 
(another piece of this kind was already subtracted during the derivation of Eq. 
([3]), see [19] for details). The first subtraction term ensures that the vacuum 
energy decreases for L — > oo. Also, we mention that this term does not enter the 
force. In this way, the renormalization done is well justified. 

We went the way from the initial mode sum (TTJ to imaginary frequencies, 
where we performed the renormalization, and back to real frequencies. 

Motivated by the initial mode sum (TTJ in Sections 2, we went the way from real 
frequencies to imaginary frequencies, where we performed the renormalization, 

3 In [22] the contributions from the photonic and from the evanescent waves were calculated. 
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and back to real frequencies. In Appendix A we demonstrated how to go the other 
way round, from imaginary to real frequencies. Here, we started directly from 
the renormalized vacuum energy (|84l) shortening the calculation significantly. It 
is needless to stress that both ways are equivalent. An interesting by-product is 
relation (11 73ft which may serve as a check for the transmission coefficient. 

As already mentioned, the two representations, one in terms of imaginary 
frequencies, the other in terms of real frequencies, coincide and that in terms of 
imaginary frequencies is just the well known one. In opposite, that in terms if 
real frequencies, has some unexpected features. First, there are the contributions 
from the waveguide modes and from the photonic modes already in the TE case. 
In the TM case, the surface modes come in addition. Second, the waveguide and 
the photonic modes show oscillations which do not decrease for large separation. 
This is clearly seen in Fig. H] for the TE case. For the TM case, the same holds 
since the difference between both cases, shown in Fig. [5J decreases for L — > oo. 

This result is somehow counterintuitive. For large separation (equally, for 
large u p ), the half-spaces become ideal conductors and the spectrum consists of 
the waveguide modes only (take, for example, u p — > oo in Fig. [2]). However, for 
any fixed L (or oj p ), we have a sum, E^+E^ t , where the oscillations compensate 
each other to a large extend. As seen in Eq. ( |89l) . the sum equals E™ ag , Eq. ( |89lh 
which deceases ~ L~ 3 . 

It must be mentioned, that the attempt to identify the contributions from 
the physical modes to the vacuum energy works well for small separation. Best 
it works in the limit L — > (non retarded case), where the surface modes of the 
TM polarization dominate the vacuum energy, 

£o~-c 2 ^f for L-»0, (154) 

(see Eqs. (11211) ). while the contributions from the waveguide and photonic modes 
of both polarizations are subleading. However, as shown in [4J, for increasing 
separation, there is a compensation between the two surface modes, and further 
increasing L, there is a compensation between the sum of the contributions from 
the surface modes and from the remaining modefl In |4], this remaining part was 
calculated as the difference between the energy calculated on the imaginary axis 
and the surface modes's contribution. This remaining part was called 'photonic'. 
But now we see that it consists of two parts, wageguide and photonic, and we 
observe also large compensations between these two. 

As a problem for further investigations we mention that the compensations 
between the waveguide and the photonic modes for large separation, possibly, 
can be handled analytically. This would make calculation of the vacuum energy 
in terms of real frequencies easier. 

In the present paper we considered the plasma model with permittivity given 
by Eq. ([2]). It should be mentioned that the separation into waveguide and 

4 The latter compensation was mentioned earlier in [3]. 
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photonic modes will appear also in the easier case of a constant permittivity (one 
needs either to take e < 1 or to consider a suitable ep in the gap). And, of course, 
this separation will show up if considering more complicated permittivities, at 
least such without dissipation. 
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A From imaginary to real frequencies 

In the main text, in Section 3, we considered the transition from real frequencies 
in the vacuum energy, to be inserted into the mode sum, Eq. (CQ), to imaginary 
frequencies in order to establish the relation with Eq. (|5]). In this appendix we 
consider the reversed transition, from imaginary frequencies as used in Eq. f!75|) . 
to real frequencies, as used in Eq. fl85|) . Of course, both are equivalent. However, 
in Section 3 we started from the unrenormalized vacuum energy, whereas here we 
can start from the already renormalized one, which makes the calculations easier. 
Further we restrict ourselves to the TE case. The corresponding formulas for the 
TM case can be written down too, but that would not bring new insights. 

We start from representation (ITS"]) of the renormalized vacuum energy in terms 
of imaginary frequencies, 



_i r°° d 

" T 

where we denoted, using ( 1761) . 



^=lW„ ^MO, (155) 



t(0 = T™(iO = — - 2 (156) 

with x = y/u)p + £ 2 . In the complex w-plane (see Fig. (jBJ)), we integrate in 
(I155p over the imaginary axis, u = i£. We note that the integrand, consid- 
ered as function of £, has poles in £ = ±iqJ E , which correspond to the fre- 
quencies P3"j) of the waveguide modes. Note, that in Eq. ()155p . the integra- 
tion over the momenta k|| is already carried out and the relation between the 
momenta and the frequency is given by Eq. ( TT3|) with formally set k\\ = 0. 
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Now we turn the integration path in 
(11551) to the right, £ — > —iu. Using £ 3 — > 
iu 3 , we get 



f->iw p 



E, 



TE 



127T 2 



x 3 d 



duu c 



du 



(157) 

lnt(-i(w +i0)), 



Figure 6: The complex cu-plane. On 
where '+ iO' indicates that the integra- the real axis, the crosses denote the 
tion path approaches the real w-axis from poles corresponding the to waveg- 
above. We separate this expression into uide modes. The dashed lines indi- 
two parts, cate ways to rotate the integration 

path. 



where 



^hnag — A + B, 



A 



12tt 2 



(158) 



duu 3 -^- lnt(-i(u + zO)) 
dui 



(159) 



has the integration region u G [0,w p ]. This region corresponds to the waves 
which are usually called evanescent. Accounting for the poles of the integrand, 
we represent this integral as the sum of the residua from passing these poles 
(entering with a factor —in since t(—iu>) has poles) and a Vp-integral, 



d 

Vp / duju! 3 —— ln£(— icu). 
' duj 



(160) 



The residua collect just into the sum over the waveguide modes, Eq. (jl9"j) . for 
s = 0. In the transmission coefficient entering this integral we note that x = 
^J^p — is real. Hence, 



can be rewritten as 



with the phase 



t(—iu) 



t(—iu) 



(*±£0 e 



-i(f> 



-2i sin(0) 



2iu)L 



1 , X + IUJ 

(p = - in h ooL. 

i x — icu 

This allows separating the Vp-integral into real and imaginary parts, 



(161) 

(162) 
(163) 



dojto 3 — 

dbJ 



Yin' 



d 

Vp / du u 3 — In (-2i sm(<j>)) . (164) 
dco 
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The integral containing the phase 0, Eq. (I163p . has no poles in its integrand and 
it can be calculated easily, 

^ 3 -r%=;rl + ^i (165) 



12tt 2 J cL; 9tt 2 48tt 2 ' 

It combines in (11641) with £7™ just into E^ E , Eq. (IH7I) . Hence, the real part of 
A is simply 

Re(A) = EH ( 166 ) 
The imaginary part of A is the remaining Vp-integral (the last term in Eq. (I164j) ) , 
which we keep as is for the moment. 

Now we turn to the remaining part of the integration region, u £ [u p , oo), i.e., 
to part B in (I158p . which we represent in the form 

I r°° d 

B = -/ duu 3 — [lm{lnt(-iu)) - iRe(lnt(-iu))} . (167) 

J Wp du 

Here we note x = — ik with k = ^Joj 2 — O0p and 

t(-iu) = T? E {k) (168) 

is just the same transmission coefficient as used in ( 1521) . In this way, we get for 
the real part of B, 

i r°° d 

Re(jB) = ~iW duJuj3 'duj lm ^ n T i ( k ) ) ' ( 169 ) 

which differs from Eq. ( 1521) merely by a change of the integration variable. So 
we have 

MB) = E™. (170) 

Collecting from (11661) and (11701) . we get from the rotation £ — > —ioj of the inte- 
gration path, 

E™ = E™ + EZ t , (171) 

which is just the same as (1851) . 

The imaginary parts, appearing in A and in B, must add up to zero, 

lm(A) + lm(B) 

= ^['^i M-« SUlW , + jj^l^ QnTH*)) • 

= 0, (172) 

since the initial expression, Eq. f l 1 5 5 j) . is real. We mention that the two integrals 
in (11721) can be joined into one using (11621) . such that the relation 

i r°° d 

dtuu 3 —Re(lnT? E (k))=0 (173) 
lztt Jo du 
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holds. 

There is also the possibility to turn the integration path the other way round, 
£ — > iu, i. e., to the left in Fig. El All calculations done above can be repeated 
with corresponding changes of signs. As a result, one obtains just the same real 
parts, i. e., Eq. ( 1171ft . The imaginary parts appear also the same as Eq. (1172j) . 
but with opposite sign. 

In this way, we demonstrated how to go the way back in the complex fre- 
quency plane as compared to Section 3. A difference is that we used in Section 
3 the variable k for the rotation, in this appendix we use the variable u, which 
appears to be more instructive. It is needless to stress that both options are 
equivalent. We obtained, as expected, the same relation (|85|) as before. The 
procedure represented in this appendix appears easier than that in Section 3. 
This is so since we consider here the vacuum energy after renormalization, Eq. 
( I83p . It is interesting to note that the subtraction terms in the renormalized 
energy of the waveguide modes, Eq. (1871) . appear also from the integral over the 
frequencies of the evanescent modes, Eq. (I165p . without explicit reference to the 
renormalization. 



B Contribution from the surface plasmons at 
small separation 

In this appendix we calculate the individual contributions from both surface 
plasmons to the vacuum energy at small separation. These dominate the energy 
and their sum is the same as given by Eq. ( 1121H . 

The contribution from a surface plasmon to the vacuum energy is given by 
Eq. (196|) with s = 5 = 0, which we consider now for u p — > 0. For this, we need 
the momenta r] s f (see Eq. (137]) ). which are solutions of Eqs. (135"1) . From this 
equation, it follows that for u p < k\\ the inequality 



'kl-ul^jKkw (174) 

holds. This motivates, for small ujp the ansatz 

T] sf = k\\ - a sf ujl + ... . (175) 
Inserting the ansatz into Eq. ( 1351) we get 



1 1 



2kn 14- nnth h. L 



1 + coth -\- 
1 



a a = 7T, hTT ( 176 ) 

2fc|| i + tanh^ 



35 



From here, using ( 139]) and (!37|) . the expansions of the frequencies follow, 

Us - OJ single = \ , 1 = ~ ~7= I + 

/l + co thM V2 

W« - ^single = | , = ~ ~7= | + O(wJ). (177) 

/l + tanh Y V 
Inserting these into fl39|) and fl37j) and integrating over k|| one comes to 

£ s = c s -^ + ... , 

E a = c a ^| + . . . , (178) 



with 



7T 2 - 24(1 + (-1 + | In 2) In 2) 
-0.030576, 

' [2(-l + V2) + (2 -In 2) In 2 



2^ 

+arcsinh(l) (2(-l + In 2) - arcsinh(l)) - Li 2 (3 - 2y/2)) , 
~ 0.0266708. (179) 

The sum of these gives Eq. (j!2ip . 
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